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■ Abstract. The paper introduces external ellipsoidal and external sphero-conal ft,- harmonics 

| for the Dunkl-Laplacian. These external /i-harmonics admit integral representations, and 

£f) • they are connected by a formula of Niven's type. External ft-harmonics in the plane are 

! expressed in terms of Jacobi polynomials P^ and Jacobi's functions Q"'^ 3 of the second 

kind. 
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1 Introduction 

In a previous article [21] the author introduced (internal) ellipsoidal /i-harmonics F n p and 
sphero-conal /i-harmonics G n>p for the Dunkl-Laplacian. The functions F n>p (x) and G niP (x) 
^f) . are polynomials in the variables x = (xq, x±, . . . , Xh) and they are /i-harmonic, that is, they 

satisfy the differential-difference equation 



A h u:= {Vl + Vl + ... + T>l)u = 0, (1.1) 
where T>j is defined by 



>< I 2> iU (x) := d 4^+ «/ (x) "^ (x)) 

■ ° x 3 x 3 

and o~j is the reflection at the jth coordinate plane. The parameters otj enter the weight function 

h(x) :=|x r|xi| ai ■■■\x k \ a K (1.2) 

Following the book [9] by Dunkl and Xu and the paper [23] by Xu we will assume that ctj > 
for each j = 0,1, ... ,k although one would expect that the range of validity can be extended 
analytically to the domain ay > — | for each j. We will also exclude the case k = 1, ceo = ol\ = 
because we want the constant fi defined below in (jl.5p to be positive. 

The parity vector p = (po,Pi, ■ ■ ■ ,Pk) has components in {0,1} and indicates that F np 
and Gn jP have parity p which means that they are sums of monomials Xq 0+2iJ0 x^ 1+2qi ■ ■ ■ x p ^ +2qk 
with qj 6 No = {0, 1,2,...}. The vector n = (n\, n^, ■ ■ ■ , n p ) counts the zeros of the corres- 
ponding Stieltjes quasi-polynomials inside k adjacent open intervals. If we set m = 2|n| + |p| := 
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2(m + • • • + rtfe) + (po + • • • + Pfe) then Gr njP is a spherical /i-harmonic of degree m, that is, it 
is a /i-harmonic homogeneous polynomial of degree m. Moreover, F n p is a polynomial of total 
degree m which, however, is usually not homogeneous. The sphero-conal /i-harmonics of fixed 
degree m form an orthogonal basis in the linear space of spherical /i-harmonics of degree m with 
respect to the inner product 



(u,v) 



/i 2 (x)u(x)f (x) dS(x), 



(1.3) 



where S k denotes the unit sphere in Equation (jl.3p uses the standard surface integral 

over the sphere so that J„ h IdS is the surface area of S fc . 

In the present paper the theory is extended to include external ellipsoidal /i-harmonics J~ n ,p 
and external sphero-conal /i-harmonics Gn,p- Throughout, we will use calligraphic letters to 
denote the "external" functions. The functions G n ,p and J~ n ,p are still solutions of (jl.ip but 
they are no longer polynomials. They have proper decay as ||x|| — ► oo which makes them 
suitable for solving problems concerning the exterior of ellipsoids or spheres, respectively. The 
external sphero-conal /i-harmonics are easy to define. They are given by 



£n, P (x) := ||x| 



-2fi-2m 



Gn,p(x), 



771 



2|n| + |p| 



(1.4) 



where 



ao + «i + h a/- + 



k-1 



(1.5) 



The definition of external ellipsoidal /i-harmonics is postponed to Section [21 

After introducing the fundamental solution of (jl.ip in Section [3] we obtain integral represen- 
tations of F^p and Gn,p in terms of the fundamental solution in Section HI Section [5] states some 
integral formulas related to Section 5.2 in [S]. As a main result of this paper, in Section we 
prove a Niven type formula that expresses ^n.p in terms of Gn,p- In the classical situation k = 2, 
a = 0, this formula was first given by Niven [TB] in an amazing award winning paper. Modified 
versions of some of Niven's result are contained in the books by Whittaker and Watson [22] and 
by Hobson [T5]. It was the attempt to extend Niven's results to external /i-harmonics that led 
to this paper. We also refer to paper [15] by Kalnins and Miller which contains a generalization 
of Niven's work from a different perspective. 

In Section [7] we illustrate the results in the relatively simple but nontrivial planar case k = 1. 
In this situation the internal ellipsoidal and sphero-conal /i-harmonics are expressed in terms 
of Jacobi polynomials pjf*'^ while the external ellipsoidal /i-harmonic ^n.p involves the second 
solution Qn"'^ of the differential equation whose first solution is P^ 1 ^. 

We should mention that there is an increasing number of papers in Applied Mathematics 
that use ellipsoidal harmonics; for example, see [3111 [5]. For applications of Dunkl operators 
and the associated Laplacian to mathematical physics, special functions, probability theory and 
geometry, we refer to [3 El H3 EH Q2] . 



2 External ellipsoidal harmonics 

We consider the Fuchsian differential equation 



k 

IK' 



k i 

a 3 -r 2 y l 



t — dj 
3=0 3 



+ 



fc-i 



j=0 



t — a 



i=0 



v = 



(2.1) 
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for the function v (t) where the prime denotes differentiation with respect to t, and 
k 



i=0 



The parameters ay are from (|1.2|) . p = (pojPi • • • jffe) is a fixed parity vector, 

ao < a\ < ■ ■ ■ < afc (2.2) 

are chosen numbers and Ao, Ai, . . . , \k~i are real spectral parameters. 

Let n = (ni, ri2, • • • , nfc) G Ng. In [2T1 Section 2] we considered the quasi-polynomial E n p 
originally introduced by Stieltjes [19] in the special case p = 0. The quasi-polynomial En,p 
has rij zeros in (aj-i, aj) for each j = 1, 2, . . . , k, and it is a solution of (|2.ip for special values of 
the parameters Ao, Ai, . . . , \k—i- F° r this set of parameters, — ?r and + ?r are the exponents of 
equation (|2.ip at infinity, where m and n are according to (|1.4p . (|1.5|) . Now ^ n ,p is the Frobenius 
solution belonging to the exponent — -y- We introduce a second Frobenius solution £ np of 
equation (|2.ip belonging to the exponent + ^ at infinity. This solution is defined for t > a/. 
and it is normalized by the condition that 

lim t^^SnJt) = 1. 

t^oo 

This normalization and the corresponding one of E n p leads to the Wronskian 

k 

SnAVKjt) ~ E n , P (-tK,p(t) = (/" + ™) ]J(t - aj)'^ (2.3) 

3=0 

which holds for t > In the classical case k = 2, a = 0, the second solution £ np appears in 
[H § 9.7] and QH § 100] . 

We introduce ellipsoidal coordinates involving the parameters (|2.2p . For every {xq, . . . ,Xk) 
in the positive cone of 

x > 0, . . . , x k > 0, (2.4) 
its ellipsoidal coordinates to, t\, . . . , tf. lie in the intervals 

dk < t < oo, ai-i<ti<ai, i = l,...,k, (2.5) 
and satisfy 

= 1 for * = 0,l,...,fc. 

J= o J 

These coordinates provide a bijective mapping between the positive cone (|2.4p and the cube (|2.5p . 

As explained in [2U Section 3], equation (jl.ip can be solved by the method of separation of 
variables in ellipsoidal coordinates. In [21] we considered the internal ellipsoidal /i-harmonic 

k 

i ? n,p(x) = n^n,p(t,), 

where (to,ii, . . . ,tk) are ellipsoidal coordinates of x. We now introduce the external ellipsoidal 
/i-harmonic ^n.p by 

^„,p(x) :=^p(<d)n = |^^n, P (x). 
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This function is originally defined in the positive cone (|2.4p and can then be extended analytically 
to the complement of the degenerate ellipsoid 

k-i x 2 

V 3 —<l, x k = 0. (2.6) 

U ak - aj 

We see this as follows. The ellipsoidal coordinate to > a k is an analytic function of x for x 
outside the degenerate ellipsoid (|2.6p . Further, E np and £ np are analytic functions on the 
interval (a k ,oo), -En.p has no zero there and -F n ,p(x) is a polynomial in x. 
We note that ^n.p has parity p, that is, 

^n,p(^(x)) = (-ir^ n ,p(x). 

Therefore, by construction, J~ n ,p satisfies equation (jl.ip . It may be appropriate here to remark 
that a solution of (jl.ip is always implicitly assumed to be defined on an open set which is 
invariant under the reflections aj, j = 0,1, ... , k. Of course, our domain of definition of ^n.p 
has this property. 



3 The fundamental solution 

We introduce the function 

k 

g(x,z):= ^ / (*(x, Z ,T))-' i nc Q .(l + r i )(l-T|p- 1 dr, (3.1) 

4 n i =or(«i + 3) [ _J ]k+1 3=0 

where 

*(x, z, t) := ||x|| 2 - 2(t x z H h T k x k z k ) + ||z|| 2 > 0, 

c u := (B(t£, with B denoting the Beta function, and fi is from (II. 5p . If aj = then delete 

the factor containing c ai , set r% = 1 in \& and omit the integration on Tj . If Q!j = for all j then 

reduces to the standard fundamental solution of Laplace's equation [lOl page 22]. 

If \xj\ \zj\ for at least one j, then \&(x, z, •) has a positive lower bound on [— l,l] fc+1 . 
Therefore, <I>(x, z) is analytic everywhere in M. k+l x ]R fc+1 except at points where \xj\ = \zj\ 
for all j. We recognize that $(x, z) = $>(z, x), 3>(7x, jz) = |7|~ 2 ^<I>(x, z) and 3>(x, z) > 0. 
Moreover, 

$(0,z) = - r(/i) ||z||- 2 ^. 

4 n T(a, + i) 

If z 7^ we can expand $(-,z) in a Taylor series at x = 0. Its coefficients are given by the 
following lemma. 

Lemma 1. If z ^ i/ien, /or ewery (go, ?i, • • • , <&) G Nq +1 , 

^•••c>S$(x,z) 



qk \ Ux Ux k 



x=0 

= Vf ■ ■ • 2?»*(x, z) | x=0 = (-1)^ • • • 2?»*(0, z), (3.2) 
w/iere |q| = q H h g fe . 
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Proof. We note that the first equality in (|3,2p is trivial so it is enough to prove the second one. 
The proof of the second equality involves some lengthy formulas, so we will use some common 
abbreviations. If j = (jo,ji, ■ ■ ■ ,jk) is a multi-index we set j! = jo! • • • j'fc!, = di° ■ ■ ■ di k k and 
x J = Xg° • • • x? k k . By Faa di Bruno's formula (or by equation |21l (5.4)]), we write 

^/(*(x, z, r)) = £ 2 l^" 2 UI/^l-|jl)(M/(x, z, r))i d^(x - • • ■ (a* - r k z k f\ 

j 

where f(u) : = and the summation is over all multi- indices j = (jo,ji, ■ ■ ■ , j k ) with j± > 
and qi — 2ji > 0, i = 0, 1, . . . , k. Computing the partial derivatives on the right-hand side and 
setting x = gives 



<2fc-2jfc 



^/(*(x,z,r))| x=0 

= (-l)N £ 2l^^l/^IHjl)(||z|| 2 ) jT ^-^(r zo) 9o - 2jo • • • (r k z k ) 

We now carry out the integrations indicated in (|3.ip taking into account that 



c aj J^(l + r j )(l-T^- 1 dr j 



Note that P™(z™) is a constant depending on j and m (the "Dunkl factorial".) We obtain 



dS&hc, z) = (-i)l q l V 2 lql- 2 ljl/(lql-ljl)(||z|| 2 ) ?! z q - 2 J, 

x=o v 7 4^ j!P^ 2j zq- 2 j 



where $ is defined by the right-hand side of (J3JJ) but with the normalization factor in front of 
the integral omitted. Now we use the first identity in (13. 2p to obtain 

PS^fx, z) 



z q-2j 



x=o ^ j! 2j z q_2 J 

j z 

= (_i)W 2 |q| " 2|j| / (|qHj|) (||z|| 2 )^ V 2 J z*. 

j J ' 

Finally, we apply |21| Lemma 2] to the function A(vq, . . . , v k ) = f(vo + • • • + v k ) and obtain the 
second identity in (|3.2p . ■ 



Since ||z||~ 2 ^ is /i-harmonic, Lemma Q] shows that $(x, •) is /i-harmonic for fixed x and so, 
by symmetry, $(-,z) is /i-harmonic for fixed z. 

Knowledge of the exact nature of the singularities of <3? will be immaterial in this paper. An 
analysis of the singularities would be along the following lines. We see that the singularity of 
$(x, z) at x = z is determined by integration over values ta close to 1 so it is convenient to 
substitute a j = 1 — Tj. If all Zj are nonzero we have the elementary inequality 

*(x,z,t) >c(||x-z|| 2 + |H|) 

if x is close enough to z and c is a positive constant; see the proof of [£l Theorem 5.5.7] for 
similar estimates. It follows that $(-,z) and its normal derivative on spheres ||x — z|| = e do not 
grow faster as e — > than in the classical case a = 0. Obvious modifications of the growth rate 
are obtained when one or several Zj vanish. If \xj\ = \zj\ for all j but Xi ^ z% for at least one i 
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then the factor 1 + n appearing in (|3.ip comes into play and makes the singularity milder than 
that at x = z. 

We now express <1> in terms of the reproducing kernel P m (x, y) of the linear space of spherical 
/i-harmonics of degree m with respect to the inner product (|1.3j) : see [9l Section 5.3]. We 
normalize the reproducing kernel so that 



y(x)= / fc 2 (x)i> m (x,y)Y(y)dS(y) 
Js k 



for every spherical /i-harmonic Y of degree m and every x G S k . The generating function for 
Gegenbauer polynomials j9j Definition 1.4.10] is given by 

oo 

(1 - 2us + s 2 )-" = C^(u)s m . (3.3) 

m=0 

The convergence is uniform for —1 < u < 1, |s| < q < 1. In (|3.3|) we set s = ||x||/||z|| and 
u = ||x|| _1 ||z|| _1 (roXo^o + ••• + t^x^z^). Then (|3.ip and the integral representation of the 
reproducing kernel [91 Theorem 5.5.5] yields the following theorem. 

Theorem 1. // ||x|| < ||z|| then 

1 ,, ,,™ ,, „ ( X Z 

The convergence is uniform for -H < q < 1. 



„_n 2 (^ + m ) Vll x ll ll z ll/ 



When fe = 2, a = 0, formula (13. 4D shows how to expand the reciprocal distance between z 
and x into a series of Legendre polynomials. This expansion is due to Laplace and it marks the 
very beginning of the history of spherical harmonics Jill page 3] . 

4 Integral formulas for external ^-harmonics 



Multiplying (|3.4p by a spherical /i-harmonic and integrating, we obtain an integral representation 
of external spherical /i-harmonics. 

Theorem 2. Let Y be a spherical h-harmonic of degree m, and let y denote the corresponding 
external spherical h-harmonic y{~x) := ||x||~ 2/i ~ 2m y(x). Then, for ||z|| > 1, 

y(z) = 2(/i + m) / /i 2 (x)$(x,z)Y(x)(iS(x). (4.1) 

In particular, Theorem [2] can be applied to a sphero-conal /i-harmonic Y = G niP . 
Let Y be any spherical /i-harmonic of degree m. If ||z|| > 1 we use Green's formula for 
/idiarmonic functions [9j Lemma 5.1.5] inside the unit ball to obtain 

0=/ ^(x)fy(x)^l- $ (x )Z )^dS(x) ) (4.2) 

where v denotes the outward normal. On the unit sphere we have 

dY dy 

—— = mY = ?ny, —— = (—2li — m)Y. 

av av 
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Therefore, adding (|4.2p to (|4.ip gives 

W=^*-(x)( W »^-.(x,.,^)« W . 

When a = this formula corresponds exactly to the well-known representation formula of 
harmonic functions in terms of a fundamental solution [10\ page 34]. By using Green's formula 
in the region between two spheres we see that we also have 

y { ,) = J fc»(x) (w^j^ - *(x, .)®^) Od) (4.3) 

when ||z|| > R, where 

S k R := {x : ||x|| = R}. 

We introduce ellipsoidal coordinates by choosing numbers ctj according to (|2.2|) . and prove 
an integral representation of external ellipsoidal /i-harmonics. 

Theorem 3. Let J denote the ellipsoid 

J {"^A '} (4 ' 4) 

for some fixed £ > a&. Then, for all z exterior to J , 

k 

^n,p(z) = yr n(t-qiP-» / / i 2 (y)^(y)^(y,z)F n , p (y)dJ(y), (4.5) 
where m = 2|n| + |p| and to is defined by 

l fe y 2 

^yT := g(t^F- (46) 

Proof. We will suppress the subscripts n, p in this proof. Both sides of equation (|4,5p are 
analytic in the exterior of J so it will be sufficient to prove the equation for ||z|| > R, where R 
is chosen so large that J fits inside the sphere S^,. 

We claim that formula (14. 3h is also true if we replace y by T . To see this we expand J- on 
the sphere in a series 

oo 

^(x) = ^y g (x), 

q=0 

where Y q is a spherical /i-harmonic of degree q. Using an uniqueness theorem for /i-harmonic 
functions defined in the exterior of spheres with proper decay at infinity (the analogue of 
Theorem 2.9]) we obtain 

oo 

^( x ) = ^i^+ 2g iixir 2 ^- 2 %(x). 

This series and the series obtained by term-by-term differentiation converge uniformly to the 
proper limits on and exterior to S^. Therefore, by taking limits we obtain 



J S k V ov* OV ) 
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Using Green's formula in the region between J and S^, we also have 

m = I k\y) (>(y)^ - *(y, *>(y), (4.7) 

where v denotes the outward unit normal vector. We now apply Green's formula to the internal 
ellipsoidal /i-harmonic F in the interior of J and obtain 

//K F(y) ^^-^ z) ^) dj(y)=o - (4 - s) 

We multiply (g2J) by E(t), (jMD by £(t) and subtract noting that E(t)J : (y) = £(t)F(y) on J. 
We obtain 



E{t)F{z) = J h 2 (y)$(y, z) (^)^^ - rfJ (y)- 



Since the ellipsoidal coordinates are orthogonal, the normal derivative is expressible in terms of 
the derivative with respect to the ellipsoidal coordinate to: 

Therefore, we find 



We replace the Wronskian by f|2.3[) and (|4.5j) is established. ■ 

In Section [6] the integral in (|4.5p will be evaluated using the formulas given in Section [5j 
Theorem [3] leads to an expansion of $ in ellipsoidal /i-harmonics. 

Theorem 4. Zei the ellipsoidal to-coordinate of z be larger than the corresponding one for y. 
Then 

*^-» = g 2 0. + ^ + |p l ) f, -*^W- < 49) 



where the positive constants e np are determined by 
e 2 ajp [ /^(x)^ ( x)dS (x) = l. 



Proof. Let J be the ellipsoid (|4.4p that contains y. The function 3>(-,z) is analytic and /i- 
harmonic on and inside J. Therefore, it can be expanded in ellipsoidal /i-harmonics as in [2H 
Section 7] and the expansion coefficients can be evaluated using (|4.5p . This gives (14, 9p after 
a simple calculation. ■ 



When k = 2, a = 0, formula (|4.9p can be found in [HI page 172]. 
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5 Integral formulas for spherical fa-harmonics 

Theorem 5. Let f be a homogeneous polynomial of degree £, and let Y be a spherical h-harmonic 
of degree m. If £ = m + 2r, r = 0, 1, 2, . . . , then 

where Y(T>) = Y(T>q, . . . ,T>k)- If £ < m or £ — m is odd then the integral in (|5.ip vanishes. 

Proof. Theorem 5.1.15 in [9] gives the expansion of / in spherical /i-harmonics on S&. It follows 
from this expansion and [U Theorem 5.1.6] that the integral in (15. ip vanishes when £ < m or 
when I — m is odd. If £ = m + 2r, r = 0, 1, 2, . . . , then [9j Theorem 5.1.15] shows that 

/ fc 2 (x)/(x)y(x)dS(x)= [ /i 2 (x)/»y(x)dS(x), (5.2) 
Js fc Js k 

where / is a homogeneous polynomial of degree m given by 

= at i? I ZTi A ^ /(x) - 
4 r r!(m + ^ + l) r 

We evaluate the integral on the right-hand side of (15, 2p by applying Theorem 5.2.4 in [9]. We 
obtain 

/ /i 2 (x)/(x)y(x)(is(x) = - — — 3— -Aiy(x>)/(x), 

y S fc c / 2 m 4 r r!(^ + l) m (m + / u+ l) r ft v >jy h 

where 

fc 

2nr( ai + i) 

/i 2 (x)dS(x) = — . 

After simplification of the constant factors we arrive at the desired equation (|5.ip . ■ 

In the special case k = 2, a = formula (|5.ip is due to Hobson [TJ1 § 100]. 
We state a consequence of Theorem [TJ 

oo 

Corollary 1. Let /(x) = ^ /^(x), where fi is a homogeneous polynomial of degree I, and the 
convergence is uniform on S k . Let Y be a spherical h-harmonic of degree m. Then 

k 

, oo nr(a, + i) 

/ /i 2 (x)/(x)y(x) dS(x) = V M0 = , '~° r A£y(X>)/ m+2r (x). (5.3) 

6 Niven's formula for external ellipsoidal fa-harmonics 

We consider the ellipsoid J from (14. 4p with semi-axes dj := ^Jt — aj . By substituting yj = djXj, 
y 6 J is transformed to x G S*. If /(x) is a continuous function on S fc , then 

*)-/(!,..,!) <a.i) 
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becomes a continuous function on J and 

k 



L 



f(x) dS(x) = ]J dj 1 f w(y)g(y) dJ(y) (6.2) 



with w(y) from (|4,6|) . 

Let Y denote a spherical /i-harmonic of degree m, and let / be as in Corollary [TJ We use the 
substitution yj = djXj and (|6.2p to write (|5.3p as 



jh\y) W {y)g{y)Y (f,---,f) dJfcO 



k oo nr(« i + i) 

= n d i 2ai+i e 2m+2 ,-i r!r ; ? ; +M+r+1) A ^(^o, . . . , ^w^, m 

j=0 r=0 

where A is the operator 

K:=dlVl + ...+dlVl 

and g is defined by (|6.ip with the same connection between and gg. We use (|6.3p in the 
special case V = G niP , m = 2|n| + |p|, and g = $(-,z), where ||z|| > do (the largest semi-axis 
of J). We also note that 

E n , P (t)G n , P (|,..-,|) =F n , p (y), 

see [2H (6.1)]. Then the integral representation (|4.5p for the external ellipsoidal /i-harmonic ^"n.p 
and ([63D yield 



0( M v oo 1 n%+j) 

2(/x + m) 1 j=o 

,pW ^n, P (i) ^^-Mrfm + ^ + r + l) 

We now use the identity 

G n , p (doV , . . . , d k V k )g(x) = E n>p (t)G n>p {T))g{x) 

which holds for any /i-harmonic function g; see the formula before (6.3) in [21]. Then we replace 
the operator —A by clqDq + • • • + a k V>\ which is possible because $(-,z) is /i-harmonic. Finally, 
we use Lemma [1] and arrive at the following theorem. 

Theorem 6. The external ellipsoidal h-harmonic J~ n)P admits the differential- difference repre- 
sentation 

^ P ( z ) = E o^Clr m) (aoP ° 2 + " " + ^'w^ii-ir^ (6.4) 

r=0 ■KH'Jm+r+l 

where m = 2|n| + |p| and \i is defined in (j 1 . 5 j) . The expansion is valid for ||z|| 2 > a k — ao- 
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Formula (|6.4p is true for all z with ||z|| 2 > a& — ao because we can choose the number t 
defining the ellipsoid (|4.4|) very close to aj.. 

We can rewrite (|6.4p slightly by using the connection between internal and external spherical 
/i-harmonics [21 1 Corollary 2]: 

This leads to the following corollary. 

Corollary 2. The external ellipsoidal h-harmonic J~n^ admits the differential- difference repre- 
sentation 

00 l — ~\) r 

^•* (z) - ^ 2^! ( ; + L + i) r ^ + • ■ ■ + ^)'s».pw * 6,5 ' 

r=0 

for ||z|| 2 > a k - ao- 

7 Examples 



We illustrate some results of this paper in the planar case k = 1. We choose ao = — 1, ai = 1 
2, /3 := ao — 5- 



and set a := ai — i, (3 := ao — \- If jPo =£>i = equation (|2.ip becomes 



(t 2 - 1) (V + |±1,' + + A , = 0. (7.1) 

This equation is satisfied by the Jacobi polynomial v(t) = P^^\f) when Ao = — n(n + a+ (3+1). 
Therefore, the Stieltjes polynomial E n $ is a multiple of P^"^\ In [21] the Stieltjes polynomials 
were normalized such that their leading coefficient equals 1. This gives 

where 

, a n n\2 n 
On = ai a ' p ' ■= 



(a + P + n + l) n 

There are similar formulas for the other three parity vectors p; compare |21[ § 2]. 

A second solution of (|7.ip with Ao = — n(n + a + j3+ 1) is Qn ] see [201 § 4.61]. Comparing 
the behavior at infinity, we see that 

f n — h C)( a '^ 

<-"n,0 — {J n^in j 

where 

, , (a,/3) 9 -n-a-/3 T(2n + a + /? + 2) 
n n T(n + a + l)r(n + /? + l)' 

Again, similar formulas hold for the other three parity vectors p. 

The planar sphero-conal coordinates r, si as defined in [211 § 4] essentially agree with polar 
coordinates xq = rcos(J), x\ = rsin0. The connection is given by cos(20) = s\. The internal 
sphero-conal /j-harmonic G n< o is equal to r 2n E nt o(si) in sphero-conal coordinates. Therefore, in 
cartesian coordinates, we obtain 



G nj0 (xo,xx) = a n (x 2 +xi) n p( a ^ 
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The external sphero-conal /i-harmonic Q n ,o is given by 

GnA^,x 1 ) = a n (xl + xiy n - a - (3 - 1 P^ ) (4^4) ■ (7-2) 

The internal ellipsoidal /i-harmonic F n< o takes the form E n ^(to)E n ^{ti) in ellipsoidal coordi- 
nates. We can express to, t\ explicitly as functions of x$, x\ as follows 

t = \ [xl + xf) 2 + (i {xl + xl) 2 - x 2 Q + xl + l) 1/2 , (7.3) 

f x = i (x 2 , + xl) 2 - ( i (xl + xl) 2 - xl + xl + l) V2 . (7.4) 

It should be mentioned that analogous formulas do not exist when k > 4 and would be very 
complicated when k = 2 or k = 3. Now internal ellipsoidal /i-harmonics can be written as 

F n , (x ,xi) = a 2 n P^\to)Pt S {ti) 

with to, t\ given in (17. 3p . (17. 4p . respectively. Moreover, external ellipsoidal /i-harmonics assume 
the form 

•T^oKzi) = a n 6 n Q( [ Q ^)(io)P^ ) (ti). (7.5) 

Formula (|7.5p is meaningful for all (xo,xi) which do not lie on the common focal line —y/2 < 
xo — V%i = of the confocal ellipses determined by a constant value of the ellipsoidal 
coordinate to- 

Formula ()6.5[) now expresses the function (|7.5I) in terms of the function (17. 2p as follows. 
Theorem 7. Let n G No- i*br every (xo,xi) S M 2 u>it/i Xq + xf > 2 we Ziaue 

00 f 1 \ r / 2 2 \ 

_ V" (V 2 - V 2 Y(r 2 + Jrn-^-ipta.jS) / X ~ X l \ f 7 fi \ 

- 2 2, r!(2n + a + p + 2)p 1^1 p o) t* + x i) U§ + 4 ) ' ( } 

where on the left-hand side to, ti are given by (|7.3p . ()7.4p . respectively. 

Obviously, we could replace T>\ — T>1 by 2T>\ on the right-hand side of (|7.6p . Formula (|7.6p 
expands the left-hand side in a series of homogeneous functions whose degree decreases with r. 
The convergence will be especially good when Xq + xf is large. The formula can be checked with 
symbolic computer algebra by replacing Xj by uxj on the left-hand side and then expanding in 
an asymptotic series at u = 00. Computer algebra confirmed the correctness of (|7.6p for small 
values of n. The natural question arises whether this formula or the more general formula (16. 5p 
can be proved more directly, that is, without using any integrals in the proof. 
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